Pulse confinement in optical fibers with random dispersion 
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Short range correlated uniform noise in the dispersion coefficient, inherent in many types of 
optical fibers, broadens and eventually destroys all initially ultra-short pulses. However, under the 
constraint that the integral of the random component of the dispersion coefficient is set to zero, or 
pinned, periodically or quasi-periodically along the fiber, the nature of the pulse propagation changes 
dramatically. For the case that randomness is added to constant positive dispersion, the pinning 
restriction significantly reduces pulse broadening. If the randomness is added to piecewise constant 
periodic dispersion, the pinning may even provide probability distributions of pulse parameters that 
are numerically indistinguishable from the statistically steady case. The pinning method can be 
used to both manufacture better fibers and upgrade existing fiber links. 



The effect of random perturbations in optical fibers 
increasingly attracts attention as the demand for the 
quality of transmission grows daily Q]. The impact 
of randomness on signal transmission in a single mode 
fiber is negative; it causes degradation of the signal and 
lowers transition capabilities j^. In particular, ampli- 
fier noise and noise in the fiber birefringence (dou- 
ble refraction) Q lead to random shifts in the pulse po- 
sition (timing jitter) and to pulse broadening, respec- 
tively. Both effects eventually cause a destruction of 
bit-patterns and eventual increase of the Bit-Error-Rate 
(BER), the most important parameter describing perfor- 
mance in fiber communications systems 

In the present paper, we consider the effect of random 
dispersion, which is, for ultrashort pulses, potentially as 
dangerous as the aforementioned effects. We, however, 
propose a way to significantly reduce the pulse deteriora- 
tion and eventually reduce the BER caused by the noise 
in dispersion by using passive (independent of pulse prop- 
erties) periodic control of the accumulated dispersion of 
the fiber link. Furthermore, the method may even pro- 
vide statistically steady propagation of the pulse along 
the fiber. 

Chromatic dispersion is an important characteristic of 
a medium and can significantly degrade the integrity of 
wave packets. In practice, chromatic dispersion is not 
uniformly distributed and often exhibits random varia- 
tions in space and time. On the other hand, wave prop- 
agation through the medium is usually much faster than 
temporal variations of the chromatic dispersion. There- 
fore, these random variations can be treated as "spa- 
tial" multiplicative noise that does not change in time. 
This multiplicative noise is conservative and the wave 
energy remains constant during propagation through the 
medium. Recently, high precision measurements of fiber 
chromatic dispersion as function of a fiber length experi- 
mentally demonstrated the significance of the dispersion 
randomness HI 



The overall chromatic dispersion in an optical fiber 
comes from two sources. The first source is the medium 
itself. The second source is the specific geometry of the 
waveguide profile. Material dispersion in the optical fiber 
is a relatively stable parameter, uniformly distributed 
along the fiber. However, waveguide dispersion is not 
nearly as stable. Existing technology does not yet pro- 
vide accurate control of the waveguide geometry of mod- 
ern fibers, where dependence of dispersion on wavelength 
is complex. As a result, the magnitudes of random varia- 
tions of fiber chromatic dispersion are typically the same 
as, or in some cases even greater than that of the mean 
dispersion [^,0. 

In the short-wavelength regime, a universal descrip- 
tion of the signal envelope in the reference frame moving 
with the packet group velocity is given by the nonlinear 
Schrodinger equation (NLS) for the complex scalar field, 
ip{z;t), see for example 

- i^^V = rf(^)atV + SV'^- (1) 

The equation is written in the dimensionless form 
Variations in the medium enter this description through 
the dispersion coefficient d{z) = ddctiz) + £,{z), which 
is decomposed into its deterministic part, rfdet(-z), and a 
random part, S.{z). Here z is the position along the fiber 
and t is the retarded time. The initial profile ■0(0; 
localized in t. We consider two different models of deter- 
ministic dispersion, both of which are standard in fiber 
optic communications. Model (A) is the case of constant 
dispersion, ddot = '^o- In the absence of noise (^(2) = 0), 
ipo{z;t) — aexp [izdo/^^] sech[t/b], where a^b^ = do (f* 
is the peak amplitude, b is the pulse width), is an ex- 
act soliton solution of (|l|). The existence of the soliton 
Q is the result of a dynamic equilibrium between dis- 
persion and nonlinearity: the two spatial scales, nonlin- 
earity zjvl = and dispersion Zd = b^/d^, coincide. 

Model (B) is the case of dispersion management (DM), 
rfdct — do ± doM [|lO|- Here dispersion is piecewise con- 
stant: positive and negative spans alternate with period 
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zdm, and < do < doM- There is no exact solution for 
the pure (no noise) Model (B), but theoretical evidence, 
confirmed by extensive numerical studies and experimen- 
tal results, indicates the existence of a breathing solution 
(DM soliton) with a nearly Gaussian shape plj-p^ . The 
localized solution here is again due to the interplay of 
dispersion and nonlinearity. In the presence of a peri- 
odic dispersion map, however, the (DM) soliton acquires 
an important characteristic, quadratic phase (chirp). In 
contrast to conventional soliton solutions, DM solitons 
can exist for zero (or even negative) values of average 
dispersion. 

Approximate scale characteristics of the dispersion 
noise present in real fibers can be extracted from ex- 
perimental results ^,|^. These results show that the 
smallest scale of noticeable change in the dispersion value 
is approximately ^ 1 — 2km. For constant dispersion 
fibers (model A), the amplifier spacing is '--^ 50 — 60km, 
and for dispersion managed fibers (model B), the pe- 
riod of a typical dispersion map is also ^ 50 — 60km. 
These scales are much longer than that of the dispersion 
variation, justifying the idealized consideration of delta- 
correlated noise for both models. Previously, the stabil- 
ity of initial pulses in the presence of the short-range- 
correlated Gaussian uniform noise Cu with zero mean, 
(Cii(-Zi)Cii(^2)) = DS (zi — Z2), was studied for both mod- 
els (A) and (B). For model (A), adiabatic theory, valid 
if the noise is weak, shows dynamical broadening of 
the pulse and its eventual destruction [|l5|. The pulse- 
broadening effects of S^u on model (B) were studied nu- 
merically and by means of a variational approach ||l5|,|l6| . 

The problem can also be addressed in the limit of 
strong noise. On short scales the nonlinearity is weak and 
propagation is essentially linear : only the phase of the 
pulse is changed by a rapidly varying dispersion, while its 
frequency spectrum is not. In this weakly nonlinear case 
the following change of variables is suggested, tp {z; t) — 

ckv exp [-Z {u;t + co^ [J^ {d{z') - do) dz'])] ^z), 
where ipu:{z) describes the pulse evolution on longer 
scales. The equation for the noise average of the slow 
filed = (^^^(z)) derived from (0) is 

(-idz + idQiu'^) ipuj = 2 / duJi^2,3H^i+^2-i^3-(^) 



X exp 



-iA J dz' (ddot - do) 




e-^'^^-/2 ^1(^2^3. (2) 



Here, A = 



and the second exponential 



on the right of (Q), which we will call the kernel, is the av- 
erage of exp [—iA ^ (z') dz'] . z^ — z, thereby setting 
the correlation scale, z^ = [A'^D]-^ ~ [u^D]-'^ ~ b^/D. 
The approximation is justified, i.e. the nonlinearity is 
weak if z^ <^ zatl- The exponential decay of the ker- 
nel with z disrupts the balance between nonlinearity and 
dispersion necessary for steady pulse propagation. 



The grim analytical conclusion is that the natural noise 
in dispersion leads to destruction of initially localized sig- 
nal in the extremes of both weak and strong noise. Nu- 
merical study of the intermediate range shows the same 
effect. 

The natural question is that of the existence of an ar- 
tificial constraint which may reduce or completely pre- 
vent this pulse broadening. We demonstrate that such 
a constraint does indeed exist, and can be readily im- 
plemented in real fibers. All that is required is that 
the accumulated dispersion, Jp dy£,{y), is set to zero, or 
pinned, either periodically or quasi-periodically with a pe- 
riod of the order or less than z^. The resulting Gaus- 
sian nonuniform noise, 5„ with zero mean is described 

by {in{y)U{z)) ^ D{5{z-y)~ T-TT^TT-) V ^nd z be- 
long to the same segment bounded by an adjacent pair 
of pinning points. Otherwise there are no correlations. 

Consider the effect of nonuniform noise, ^ — 
^n, in the weakly nonlinear case. in (^ is 

z {Ij+i — z) / {Ij+i — Ij) with z € [Ij; Ij+i]. Thus the de- 
cay of the kernel in (^) is replaced by oscillatory (or 
quasi-oscillatory, when the period is fluctuating) behav- 
ior. In the nonuniform case, additional averaging over 
a single pinning leg reduces the z-dependent kernel of 
(I) to exp [~A^Dl/4] ^tt/ [A2L»;]Erfi [VA2i^/4j . In 

Model (B), independent averaging of the first exponential 
on the right-hand side of (||) over zdm (see also []TT| , p^ ) 
replaces it by 2sin [AzDAf/4] / [Az^Af]- The overall re- 
sult of the averaging procedure in the case of ^ = ^„ is 
the emergence of a nonlinear term which does not vanish 
as z — > 00, in contrast to the case ^ = ^u, where the 
nonlinearity dies away with z 00. One also concludes 
that zj sets the critical pinning period: the pinning is 
efficient only ii I — li^i — li z^. It is possible to show 
that the averaged equation for ^ = ^„ has a steady so- 
lution Notice also that the averaged kernel decays 



exponentially as A 



This is much faster than the 



algebraic decay obtained from averaging over the disper- 
sion map period in pure dispersion management. There- 
fore, the model of [ p^ corresponding to a very narrow 
type of nonlinear kernel (in cj-space) is better suited for 
the case considered here than for the case of pure DM for 
which it was originally proposed. The averaged equation 
approach is d priori applicable for a large but finite z, so 
the emergence of a steady localized solution for the aver- 
aged equation does not mean that the original problem 
possesses a steady state. In the direct numerical simu- 
lation of (|l|), to which we now switch our attention, the 
effect is seen as essentially limiting the pulse broadening. 

We perform numerical investigations of both models 
(A) and (B) with ^ — ^„ and ^ = ^„ in intermediate 
case, zjvL ^ z^. Fourier split-step scheme with 2^^ tem- 
poral Fourier modes and periodic conditions imposed on 
the boundaries of the domain t e [—180, 180] is imple- 
mented. The spatial step is Zstep — 0.01 and the numer- 
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ical convergence was checked by varying the size of the 
periodic box and number of the Fourier harmonics. Pa- 
rameters for the initial signal were chosen to be dp = l, 
a = 1 in model (A), and do = 0.15, doM = 0.1, zdm = 1, 
\^p{0;t)\ = 0.79exp(-tV2.6) in model (B). The setup in 
model (B) is borrowed from |Q and corresponds to ex- 
perimentally available DM fibers. Gaussian zero-mean 
noise correlated at z„ = 0.1 with amplitude dn — 1 mod- 
els the ^-correlated uniform noise with D = d^Zn — 0.1. 
The nonuniform noise is constructed from the uniform 
noise by the following subtraction at every pinning leg: 
^niz) = £,u{z)-j--^^ fi'.^^ dyi{y). Pinning strategies of 
two types are considered: strictly periodic, Ij+i — Ij = I, 
where / is fixed; and quasi-periodic, Ij^i — Ij — 1(1 + 77), 
where 77 is a random number uniformly distributed be- 
tween ±1/2. The averaged (or otherwise strict) pinning 
period for the nonuniform case is taken to be 1, 5 or 10. 
The simulation runs until the pulse arrives at z = 95. 
Statistics were collected for 10^ — 10^, and in a special 
case 10^, realizations. 



Model A, D=0.1, N=1000 
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The averaged pulse- width (full width at half maximum 
amplitude) as a function of z is shown in Figure 1, see 
also ijl^ . The capital letter subscript of the figures corre- 
sponds to the type of model, (A) or (B). Solid black, red, 
green, and blue represent uniform noise, and nonuniform 
noise with pinning period / = 1, 5, 10 respectively. The 
quasi-periodic curves are dashed and of the same color 



as the respective periodic ones. 

For Model(A), all types of nonuniform noise demon- 
strate a significant reduction in the rate of pulse broad- 
ening when compared with the uniform case. The indi- 
vidual configurations that degrade (through pulse split- 
ting, etc.) in the uniform case maintain pulse integrity 
when each type of nonuniform compensation is applied. 
The dependence on the pinning period is monotonic: the 
peak amplitude of the pulse decays faster as the pinning 
period increases. The difference between periodic and 
respective quasi-periodic cases is minor, with a slightly 
better confinement observed for the quasi-periodic case. 
The destruction of the pulse is also accompanied by emis- 
sion of continuous radiation by the soliton. The radiation 
is clearly seen in the movie made for individual runs . 
Once the radiation reaches the boundaries of the box, it 
reflects and starts to interfere with the still localized solu- 
tion. The latter shows up in the change of the averaged- 
width behavior at larger distances, z ~ 20. In principle, 
one could introduce absorbing boundary conditions to 
enable longer-time simulations. However, these condi- 
tions are artificial, as in real data streams pulses are not 
isolated. Each pulse emits continuous radiation which 
eventually interacts with neighboring pulses. In this pa- 
per we purposefully study self-interaction of a pulse and 
its radiation via refiecting boundary conditions, which is 
a simple model for the behavior of a real bit stream. Our 
results reveal that after z ~ 20, the rate of pulse disinte- 
gration is reduced, suggesting that in a real bit pattern, 
pulse disintegration might be prevented by radiation ab- 
sorption from neighboring pulses. 

The effect of nonuniform noise is more dramatic in the 
case of Model(B). For nonuniform compensation with the 
averaged period / = 1 (and also less) one observes a ten- 
dency toward statistically steady behavior: the average 
pulse width does not decay (in contrast to a decay in 
the uniform case), and the PDF of the pulse width (and 
of other variables characterizing the pulse propagation, 
such as amplitude) does not change shape with z 
There is no visible emission of radiation by the localized 
solution for any case of Model (B). We have also checked 
that temporal and spatial averages (for example for the 
PDF of the pulse width) coincide. The dependence on the 
type of compensation (for Z > 1) is monotonic, and the 
difference between the random and quasi-random cases 
is minor, as in the case of Model(A). 

Notice that the tremendous reduction in the pulse de- 
cay is achieved by very minor changes in dispersion (see 
jigf , for the comparable plots of the dispersion profiles). 

Our analysis gives practical recommendations for im- 
proving fiber system performance that is limited by ran- 
domness in chromatic dispersion. The limitation origi- 
nates from accumulation of the integral dispersion. The 
distance between naturally occurring nearest zeros of 
the accumulated dispersion grows with the fiber length 
(~ yfz as it is shown in ||l^ ) . This growth causes asymp- 
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totic decay of the nonlinearity with z and, therefore, 
pulse degradation. We have shown that the signal can be 
stabilized by periodic or quasi-periodic pinning of the ac- 
cumulated random dispersion pl| ] . This can be achieved 
by first measuring the mismatch between the nominal 
dispersion of the fiber line and the actual accumulated 
dispersion, and second, inserting a small fiber to com- 
pensate this mismatch [ p2[ . 
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